PREPRINT (December 13, 2011) 



Coexistence of Superconductivity and Magnetism in FeSei_a; under Pressure 

M. Bendele,i'2'[2] A. Ichsanow/ Yu. Pashkevich,^ L. Keller,"* Th. Strassle,"^ 
A. Gusev,'^ E. Pomjakushina,^ K. Conder,^ R. Khasanov,^ and H. Keller* 

^ Physik-Institut der Umversitat Zurich, Winterthurerstrasse 190, CH-8057 Zurich, Switzerland 
^Laboratory for Muon Spin Spectroscopy, Paul Scherrer Institut, CH-5232 Villigen PSI, Switzerland 
''A. A. Galkin Donetsk Phystech NASU, 83114 Donetsk, Ukraine 
'^Laboratory for Neutron Scattering, Paul Scherrer Institut, CH-5232 Villigen PSI, Switzerland 
^Laboratory for Developments and Methods, Paul Scherrer Institut, CH-5232 VilUgen PSI, Switzerland 

An extended investigation of the electronic phase diagram of FeSei-a; up to pressures of p ~ 
2.4 GPa by means of ac and dc magnetization, zero field muon spin rotation (ZF /xSR), and neutron 
diffraction is presented. ZF /iSR indicates that at pressures p > 0.8 GPa static magnetic order occurs 
in FeSei-a; and occupies the full sample volume for p > 1.2 GPa. ac magnetization measurements 
reveal that the superconducting volume fraction stays close to 100% up to the highest pressure 
investigated. In addition, above p > 1.2 GPa both the superconducting transition temperature Tc 
and the magnetic ordering temperature Tn increase simultaneously, and both superconductivity 
and magnetism are stabilized with increasing pressure. Calculations indicate only one possible 
muon stopping site in FeSei-a,, located on the line connecting the Se atoms along the c-direction. 
Different magnetic structures are proposed and checked by combining the muon stopping calculations 
with a symmetry analysis, leading to a similar structure as in the LaFeAsO family of Fe-based 
superconductors. Furthermore, it is shown that the magnetic moment is pressure dependent and 
with a rather small value of /i ~ 0.2 /ib at p ~ 2.4 GPa. 



PACS numbers: 76.75.+i 74.25.Dw 74.62.Fj 74.70.Xa 

I. INTRODUCTION 

Shortly after the discovery of superconductivity in the 
Fe-based compound LaFeAs0i_2;Fj^ in 2008 by Kami- 
hara et al.^ Hsu et al^ observed superconductivity in 
the basic binary compound FeSci-^:. This simple sys- 
tem shares the superconducting layers consisting of a Fe 
square planar sheet tetrahedrally coordinated by As/P 
or Se/Te atoms as a common feature with all of the Fe- 
based superconductors. Most of the known Fe-based su- 
perconductors are made up of a stack of the electron- 
ically active layers, separated by layers that act as a 
charge reservoir to dope the Fe-As/Se layers. FeSei-^^ 
is an exception to that rule because it consists of a stack 
of superconducting layers only. In this binary system 
the superconducting transition temperature is Tc ~ 8 K. 
Thus, it could be argued that this is more a conventional 
superconductor.!^ Shortly after, the electronic and mag- 
netic phase diagram under pressure was studied.ISEllt was 
found that the transition temperature exhibits one of the 
largest pressure effects on known. It reaches values of 
Tc « 37K at p « 9 GPa, demonstrating that FeSei_:r 
in fact is a high temperature superconductor. Further- 
more, it was found that tetragonal FeSci^x undergoes a 
structural phase transition starting at p ^ 9 GPa from 
a tetragonal to a hexagonal, non-superconducting and 
more densely packed phase. With increasing pressure 
the volume fraction of the tetragonal phase as well as Tc 
decrease until at high pressures {p > 20 GPa) only the 
non-superconducting hexagonal phase is present.'^ Early 
muon spin rotation (/iSR) experiments on FeSei_j; re- 
vealed that the system is non-magnetic at ambient pres- 
sure down to T = 0.02 K.l^ The investigation of the 



pressure dependence also did not show magnetic order 
in the beginning up to the highest pressures, just be- 
fore the structural phase transition occurs.l^ This is in 
striking contrast to the other Fe-based superconductors 
that usually exhibit static magnetic order in the parent 
compound. This is unexpected, since the FeSei_a; lay- 
ers are isoelectric to those of the parent compounds of 
other Fe-based superconductors.'^ Shortly after, however, 
NMR studies showed a wipeout of the signal that revealed 
an incipient magnetic phase transition under pressure.!^ 
This possibly may be interpreted as static magnetic or- 
der with a broad field distribution or as slow spin fluc- 
tuations, since no magnetic order was observed by fiSK 
at ambient pressure. It seems that both the magnetic 
and the superconducting states stabilize with increasing 
pressure. In fact, static magnetic ordering was observed 
above p ~ 0.8 GPa by means of /xSR.I^ The experiments 
revealed that as soon as magnetic ordering occurs, the 
magnetic and the superconducting states seem to com- 
pete with each other. This is because the incommensu- 
rate magnetic order gets suppressed when superconduc- 
tivity sets in and, in addition, Tc decreases in the pres- 
sure region 0.8 < p < 1.2 GPa. Above p ~ 1.2 GPa both 
ground states apparently coexist on an atomic length 
scale. Both the magnetic ordering temperature Tn and 
Tc increase simultaneously with increasing pressure, and 
the magnetic order becomes commensurate.* 

In this paper an extended study of the electronic and 
magnetic properties of FeSei_a: under pressure investi- 
gated by means of ac susceptibility and fiSK is presented. 
In addition, magnetic structures of FeSci-^: under pres- 
sure are proposed and checked by neutron diffraction 
measurements. The magnetic moment in the ordered 
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state is estimated for different pressures. Furthermore, 
the discrepancy between Mossbauei'^ and iiSR resultJ^ is 
discussed under the aspect that the samples used i n ea ch 
study were prepared by slightly different methodsPiSl 



II. SAMPLES 

The FeSei_2; samples were prepared following the pro- 
cedures described in Refs.PandTUI In both methods the 
samples are placed in sealed silica tubes and are prepared 
in two steps. In the first step Pomjakushina et al^ used 
selenium and iron powders as starting materials and syn- 
thesized FeSci-a; in a solid state reaction at temperatures 
ranging from 400 — 700° C. After powderizing the sam- 
ples in He-atmosphere, they were reannealed at 700 °C, 
then the temperature was stabilized at 420 °C, and finally 
they were cooled slowly to room temperature. McQueen 
et al.^ on the other hand, used shots of selenium iron 
pieces. They were molten at 1075 °C, powderized and 
annealed again at T ~ 400 °C. However, the main differ- 
ence of the two procedures is that the samples prepared 
by the method of McQueen et al^ are quenched from 
^ 400 °C to ~ — 15°C, whereas the samples prepared af- 
ter Pomjakushina et a^l^are cooled slowly from ^ 400 °C 
to room temperature. Here, the specimens are denoted 
as FeSeo.98 for the slowly cooled ones, and '^FeSeo.gs for 
the quenched ones. All samples were found to be phase 
pure with a superconducting transition temperature of 
Tc ~ 8K. In fact, the transitions to the superconducting 
state is for both preparation procedures very similar (see 

III. SUPERCONDUCTING PROPERTIES 

The superconducting properties of FeSei_a; were stud- 
ied by means of ac and dc magnetization measurements 
(Fig. [l]). The zero field cooled dc measurements, pre- 
formed in a commercial Quantum Design MPMS SQUID 
7T magnetomenter m iiqH — 0.2 mT, revealed T^. ~ 8K 
for both FeSeo.98 and '^FeSco.gs. The ac magnetiza- 
tion measurements under pressure were performed in a 
home made ac susceptometer in piston-cylinder pressure 
cells, especially designed for fiSH experiments. The ac 
amplitude was /xo-ffac ~ 0.1 mT and the frequency was 
I'ac = 94 Hz. As a pressure transmitting medium 7373 
Daphne oil was used. The pressure applied was measured 
in situ by monitoring the the shift of T^. of Pb or/and In. 
To ensure that the position of the sample in the cell is the 
same for all pressures investigated the pick up and exci- 
tation coils were directly wound on the pressure cell. Ad- 
ditional ac magnetization measurements were performed 
to check whether the ac signal under pressure was en- 
tirely determined by the bulk Meissner response of each 
grain. Thus, other effects like e.g. weak links between 
the individual grains or surface superconductivity can be 
excluded. This was done on a commercial Quantum De- 
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FIG. 1: (color online) (a) Temperature dependence of the nor- 
malized magnetization M/M(2K) for FeSeo.98 and '^FeSeo.gg 
(see text for details on sample notation) at p = GPa. The 
transition temperature Tc of both samples, obtained by the 
intersection of straight lines fit to the data above and be- 
low the transition is ~ 8 K. The shapes of the magnetization 
curves for the two samples are very similar, (b) Dependence 
of the ac magnetization Mac on the ac field amplitude /io-ffac 
at a fixed frequency v^c = 100 Hz (c) and the ac frequency //ac 
at a fixed ac field amplitude /io^^ac = 0.1 mT. 



sign PPMS in various fields (0 > ^iqHj^q > 0.5 mT) and 
frequencies {0 > > 599 Hz). As shown in Fig. [Tjo and c 
the experiments reveal that the ac magnetization scales 
linearly with the field and is independent of frequency as 
expected for a superconductor in the Meissner state. 

The superconducting transition temperature of 
FeSci-x (FeSeo.98 and '^FeSeo.gs) is Tc ~ 8K at ambient 
pressure (see Fig. [T]). Upon applying hydrostatic pres- 
sure FeSei-i: exhibits one of the highest pressure effects 
known on . The overall increase of is non monotonic 
and shows a local maximum at p ~ 0.8 GPa, followed 
by a local minimum at p ~ 1.2 GPa (Fig. [2^). This 
behavior is similar to that already observed earlier both 
by dc and ac magnetization.'^E^^ In the region where 
decreases static magnetism develops in the sample and 
competes with superconductivity (see below and Ref . [8| . 
Upon increasing the pressure above p = 1.2 GPa the 
superconducting transition temperature increases again 
and reaches values of ^ 16K at the highest pressure 
investigated in this study (2.4 GPa). 

In Fig. [2]d the diamagnetic response at T = 6 K nor- 
malized to the value at ambient pressure is shown as a 
function of pressure. At ambient pressure the value of the 
measured ac voltage of the samples in the pressure cell 
is equal to the magnetization of the sample measured in 
the PPMS magnetomenter without a pressure cell. Thus, 
the ac voltage in the pressure cell is representing the su- 
perconducting response of FeSei_a;. Calculating the sus- 
ceptibility from the magnetization measurements allows 
to estimate the superconducting volume fraction. The 
susceptibility was determined to Xac — 1-3 (Fig. [T^). By 
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FIG. 2; (color online) (a) Dependence of the superconduct- 
ing transition temperature on pressure p of FeSei-^;. The 
line is a guide to the eyes, (b) Pressure p dependence of the 
ac susceptibility x normalized to the ambient pressure value 
x(OGPa) at T = 6K, indicating bulk superconductivity for 
all pressures investigated. See text for details on sample no- 
tation. 



assuming the samples consist of individual sphere like 
shaped grains with a demagnetization factor of n ~ 1/3 
leads to an ideal diamagnetic response of x = — 1. This 
indicates that FeSei_a; is a bulk superconductor with a 
superconducting volume fraction close to 100%. Since 
the absolute value of the ac response measured at 6 K for 
each individual pressure is similar, it is concluded that 
the sample is a bulk superconductor up to the highest 
pressure investigated. 



IV. MAGNETIC PROPERTIES 

The magnetic response of FeSei_j; for various pres- 
sures was studied by means of zero-field muon spin rota- 
tion experiments (ZF /iSR). The experiments were car- 
ried out using the //El beam Hne at the GPD instru- 
ment at the Paul Scherrer Institute (PSI, Switzerland) at 
temperatures ranging from 0.25 to 80 K. The /iSR time 
spectra were analyzed using the free software package 
MUSRFIT.i'* ZF /iSR is a well known technique to study 



magnetically ordered phases where the muon acts as a lo- 
cal magnetic microprobe. Positively charged muons are 
implanted into the sample where they thermalize after a 
short time (< 10"^'^ s). Once stopped at an interstitial 
site the muon interacts with its local environment and 
decays after its lifetime of = 2.197/is into a positron 
and two neutrinos. The positron is emitted predomi- 
nantly along the muon spin direction at the time of de- 
cay. Thus, by monitoring the time evolution of the muon 
spin polarization, information on the local magnetic field 
at the muon stopping site i?i„t and the magnetic volume 
fraction are obtained. 

The /iSR signal in a pressure cell consists of a super- 
position of two components, one arising from the sample 
{^) and one from the pressure cell (^^'^): 



(1) 



In the data analysis the ratio of the component of the 
pressure cell and the component of the sample J^'^ j 
was kept constant for each individual pressure and was 
always w 50%. For the present study two different pres- 
sure cells consisting of MP35N and CuBe were used. The 
ZF response of the empty cells is described elsewhere.!^ 
As we reported earlier,!^ in the low pressure region, 
where increases linearly with p, no magnetic order is 
observed in all of the samples. The /xSR time spectra 
are overlapping for all temperatures, indicating the same 
magnetic state for all temperatures measured (Fig. [3^ 
and b). The /.iSR time spectra were analyzed using a 
single exponential decay function: 



JSc^{t) = ylgexp[-Aot] 



(2) 



Here Aq is the Lorentzian depolarization rate. The expo- 
nential behavior at low pressures indicates the presence of 
diluted and randomly distributed and oriented magnetic 
moments in the sample volume which can be attributed 
to traces of Fe impurities.^ 

As shown in Fig.jsj: and d for p > 0.8 GPa spontaneous 
muon-spin precession is observed, reflecting the appear- 
ance of static magnetic order below the Neel temperature 
Tn- The analysis was made by taking into account that 
the magnetic order appears gradually: one part of the 
muons experiences a static local field and the other part 
stops in non-magnetic regions: 



I 1 

-Jose exp[-Att] + - exp[-Ait] 



+ (1 - m) exp[-Aoi] 



(3) 



Here m is the magnetic volume fraction of the sample, 
/osc represents the magnetic signal of the sample and has, 
depending on pressure, the form /osc = cos(woi + </'o) 
or /osc = joi^ot + 4>o), whereas ujq is the precession fre- 
quency, Jo is a zeroth-order spherical Bessel function, and 
00 the initial phase of the muon ensemble. The param- 
eters At and Aj describe the relaxation transverse and 
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FIG. 3: (color online) . Zero-field /xSR time spectra of FeSeo.gs 
for (a) p = and (c) p = 1.4 GPa, and "^FeSei-^ for (b) p = 
and (d) p = 1.4 GPa for different temperatures. The lines are 
fits of Eq. to the data. 



longitudinal to the muon spin of the magnetic signal, re- 
spectively. 

In the pressure region where Tc decreases, both the 
magnetic and superconducting ground state are compet- 
ing. This is seen first by the decrease of Tc (Fig. [2^) 
and second by a decrease of the frequency and the mag- 
netic volume fraction m below T^. (see Fig.|4^ and b, and 
Ref. IH])- In this region the magnetic signal is described 
best by a Bessel function which indicates the presence of 
incommensurate magnetic order in the samples.'^ 

As shown in Fig. [4] for p > 1.2 GPa (where the local 
minimum of Tc is reached) superconductivity and mag- 
netic order coexist in the full sample volume. Here, the 
magnetic volume fraction reaches 100% and stays con- 
stant in the superconducting state down to the lowest 
temperature where also the superconducting volume frac- 
tion remains constant at ~ 100% (see Fig.[2|. Moreover, 
Bint is not significantly changing (decreasing) below Tc, 
and the magnetic order changes from an incommensurate 
to a commensurate as reflected in the /iSR line shape 
which is described better by a damped cosine function 



with zero initial phase than by a Bessel function. This in- 
dicates coexistance of superconductivity and magnetism 
in the full sample voulume. 

To determine the zero-temperature value of i3int(0) 
and Tn the temperature dependence of Biat (T) was fitted 
to the power law expression: 



A„t(T) = Bi„t(0) (l - ( ^) 



as P 



(4) 



Here a and /3 are the power exponents. For the pressure 
region in which Bint decreases in the superconducting 
state, only the data above Tc were used to analyze the 
data with Eq. ([4|. The obtained values of i?int(0) and 
Tn are plotted in Fig. [5^ and b together with the results 
from earlier studies of FeSeo.94 and "^FeSeo.gs!^^ For all 
samples i3int(0) increases with increasing pressure (see 
Fig. [5^) . As shown in Fig. [sj^ the Neel temperature in- 
creases in parallel from Tn = 17K atp = 0.8 GPa where 
magnetism appears in FeSci-a, with increasing pressure 
to Tn = 55K at the maximum pressure p ~ 2.4 GPa in- 
vestigated here. No tendency for a saturation at high 
pressures of both i?int(0) and Tn is observed. 

Unlike the /j,SR experiments presented here, an ear- 
lier Mossbauer study did not reveal magnetic order un- 
der pressure in FeSei-^;.^ However, the samples used in 
this study were prepared after the method proposed by 
McQueen et al!^ As mentioned already above, samples 
denoted as "^FeSei-^ were prepared following exactly 
the recipe of McQueen et alP^ and were investigated 
by means of /iSR.^ In contrast to the earlier study of 
Ref. [in] they also show a similar magnetic behavior as 
the samples prepared by our method (see Fig.|3j3 and d). 
In particular they also show magnetic order upon apply- 
ing pressure. A simple explanation of this discrepancy 
could be that magnetism was overseen. Low tempera- 
ture (T = 4.2 K) Mossbauer spectra were taken only at 
few pressures: At ambient pressure no magnetic order in 
agreement with the fj,SR experiments was seen, and at 
p = 14.4 GPa and 19.7 GPa no magnetic hyperfine split- 
ting in the Mossbauer spectra was observed. 

In Fig.jsj; i3int(0) vs. Tn is plotted, indicating that the 
magnetic moment is increasing with increasing Tn . This 
points to a more robust magnetic order with increas- 
ing pressure. When the magnetic order is fully estab- 
lished (above p = 1.2 GPa; the magnetic volume fraction 
reaches 100%) Tc starts to increase again (Fig. [2^) simul- 
taneously with Tn up to Tn ~ 60 K and Tc ~ 16 K at the 
highest investigated pressure in this study. It seems that 
both order parameters are stabilized at high pressures: 
(i) Both Tc and Tn increase with increasing pressure, (ii) 
the magnetic and superconducting volume fractions stay 
100 % even below Tc to the highest investigated pressure, 
and (iii) the internal magnetic field i?int (0) increases with 
increasing pressure for all samples studied. 
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FIG. 4: (color online) Temperature dependence of (a) the 
internal magnetic field at the muon stopping site -Bint and 
of (b) the magnetic volume fraction for FeSei_a; for various 
pressures. Both parameters are obtained directly from the fit 
of Eq. (|l| to the data. The solid lines in (a) correspond to fits 
of Bint(T) in the region Tc > T > Tn to Eq. Q. For details 
see text. The solid lines in (b) are a guide to the eyes. The 
arrows indicate the superconducting transition temperature 
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FIG. 5: (color online) (a) Pressure dependence of the internal 
magnetic field at the muon stopping site i3int(0). (b) Pres- 
sure dependence of the magnetic ordering temperature Tn. 
The dotted lines in (a) and (b) are guides to the eye. The 
dependence of Bint on Tn is shown (c). See text for details 
on sample notation. 



V. MUON STOPPING SITE AND MAGNETIC 
MOMENT 

Up to now it is not clear what kind of magnetic struc- 
ture develops in FeSei_a; under pressure. Calculations of 
the muon stopping sites at different pressures were per- 
formed and combined with a symmetry analysis to check 
for possible different magnetic structures. 

The space group symmetry of FeSei_j: at low temper- 
atures is Cmma with Fe in the 4a-position (1/4, 0, 0) and 
Se in the 4(7-position (0, 1/4, z) (see for instance Ref . IT71) . 
Here the symmetry of the FeSei_a: layers exactly resem- 
bles the symmetry of the FeAs-layers in the LaFeAsO 
compound with the same Cmma space group which re^ 
mains unchanged in FeSei_a; up to a pressure p « 9 GPa.** 

In order to evaluate possible muon sites the modified 
Thomas Fermi approach^^ and available structural data 
were used.'^ This method allows to determine directly the 
self consistent distribution of the valent electron density 



from which the electrostatic potential is obtained. Lo- 
cal interstitial minima of this potential serve as stopping 
sites for muons. The applicability of this approach was 
verified by comparing the numerical results with the ex- 
perimentally determined muon sites in EFqO^^ (R — 
rare earth) and by a successful interpretation of /iSR 
spectra of the complex magnetic structures in layered 
cobahites i?BaCo2 05.^ and Fe-pnictides iJFeAsO.I^ 

Only one possible muon stopping site is observed. It is 
located on the line connecting the Se - Se ions along the 
c-direction with the coordinates (0, 1/4, z) and has the 4f; 
local point symmetry (mm2) i.e. the same as the Se ions. 
The position of the muon sites in the crystallographic 
cell is shown in Fig. |6] Note that the crystallographic 
unit cell differs from the primitive cell which is built 
by primitive translations a\ — {a/2,b/2,Q) — {T^,Ty,0), 
02 = (-a/2,6/2,0) = {-T,,Ty,0), and 03 = (0,0, c) = 
(0,0, 2t,). 

As seen in Table |Tj application of pressure leads to a 
general increase of the distance of the calculated muon 
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FIG. 6: (color online) The crystallographic unit cell of 
FeSei-i in the Cmma setting. The enumeration of the Fe 
atoms and the muon positions is shown. 



TABLE I: The pressure dependence of the calculated muon 
position and Fe-Se-Fe bond angles aa (along the a-direction) 
and Ob (along the 6-direction) . The crystallographic data are 
from Refs.l4land[rfl 
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stopping sites to the iron afo-plane, whereas the angles 
of the Fe-Se-Fe bonds aa (along the a-direction) and ab 
(along the 6-direction) are almost identical at ambient 
pressure. However, at higher pressures they tend to dif- 
fer. 

The stronger reduction of the c-axis compared to the 
a- and &-axis leads to an increase of the Fe-Se-Fe bond 
angle that can be interpreted as a tendency to antifcrro- 
magnetic exchange in accordance with the semi empiri- 
cal Goodenough Kanamori rules.'^^HSll Note that already 
small variations of the Fe-As-Fe bond angles along a- 
and 6-axes in the i?FeAsO compounds lead to a dras- 
tic change of the magnetic exchange sign from anti- 
ferromagnetic (positive) along a-axis to ferromagnetic 
(negative).-- However, opposite to the i?FeAsO the b- 
axis remains in FeSei-^^ larger than the a-axis for all 
pressures. Due to this similarity one can suppose the oc- 
currence of a ferromagnetic type of order along the a-axis 
and an antiferromagnetic one along the 6-axis in FeSci-^: 
under pressure. The minimal model which could account 



for this feature should include a doubling of the primitive 
cell along the 6-axis with magnetic propagation vectors 
either Kj = (0, Tr/xy, 71/2x3) or Kjj = (0,7r/Ty,0). Ad- 
ditionally, more simple possible magnetic vectors such as 
Kq = (0,0,0) and Kjjj — (0, 0, 7r/2T2) are considered. 

The calculations of the symmetry analysis and the 
magnitude and symmetry of the dipole fields of the Fe 
subsystem at the muon are more rigorously discussed in 
the Appendix. Application of pressure leads to an in- 
crease of the magnetic field at the muon stopping site 
as observed in the experiments (see Fig. [s]) only for 
the Kj and Kjj translation symmetries. For the Kq 
and Kjjj translation symmetries application of pressure 
would lead to a decrease of the magnetic field. This be- 
havior can be explained as the result of a competition 
between a general constraint of the lattice constants and 
a simultaneous shifting of the muon positions further 
away from the Fe a6-plane. As a result from this fea- 
ture and the above mentioned similarity to the i?FeAsO 
family it may be concluded that only the Kj and Ku 
translation symmetries are possible symmetries of the 
magnetic structures for FeSei_j; under pressure. Com- 
paring both possible magnetic structures Ki and Kji 
(shown in Fig. [7| with the experimental data presented 
in Fig.[5]leads to magnetic fields along the z-coordinate of 
B^{Ki) = 354:.6-my{Ki) and B^{Kii) = 334:.3-my{Kii), 
respectively. Here rrii are the iron magnetic order pa- 
rameters (see Eq. ([s]). This corresponds to a Fe mag- 
netic moment /i ~ 0.2/iB for both magnetic structures. 
However, the very modest shift of the muon position in 
the region 0.25 — 4 GPa calculated here, cannot explain 
the giant increase (four times) of the internal magnetic 
field Bint with an increase of the pressure from 1 GPa to 
2.4 GPa. Therefore, all these changes are connected with 
a pressure induced increase of the iron magnetic moment. 
The right scale of Fig. [Sj; shows the estimated value of 
the magnetic moment using dipole-dipole calculations for 
4 GPa. 



VI. NEUTRON DIFFRACTION 

Neutron diffraction experiments were performed on 
the Cold Neutron Powder Diffractometer DMC at SINQ 
(PSI) at a pressure of p = 4.4(5) GPa in a Paris- 
Edinburgh presJ^ in order to investigate the proposed 
magnetic structures of FeSei_j^ on polycristalline samples 
of 40 mm'^ effective volume in the beam. The pressure 
was determined by the known pressure dependence of 
the c-axis of FeSei-^.'^ The experiments were performed 
at temperatures of 5 K and 150 K using neutrons with 
a wavelength of A = 2.4575 A. The FULLPROF pro- 
gram was applied to analyze and to model the diffraction 
data.l23 

The diffraction patterns measured at T — 5 K and 
150 K were normalized to each other, and then subtracted 
from each other in order to obtain evidence of possible 
magnetic Bragg peak. However, no difference peak was 
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FIG. 7: (color online) Possible magnetic structures of FeSei_a; 
under pressure: (a) my{Ki)-type and (b) my{Kij)-type. mi 
are the the iron magnetic order parameters (see Eq. ([5|. 



observed, except at the positions of the nuclear peaks 
(see Fig. [8^). The different intensities of the nuclear 
peaks at the investigated temperatures result from the 
temperature dependent Debye- Waller factors. There are 
two possible explanations that no magnetic Bragg peaks 
were observed with neutrons in contrast to /iSR which 
shows static magnetism: (i) the magnetic moment is too 
small, resulting in an intensity of the magnetic diffraction 
peak that is hidden below the background of the sample 
and the pressure cell, or (ii) the magnetic order is static, 
but no long range order occurs (muons are sensitive only 
over a few unit cells). 

However, because oscillations are seen in the /iSR time 
spectra (see Fig. [s]) the magnetic order is long range, 
thus leading to the conclusion that static magnetic order 
occurs below Tn. The muon stopping site calculations 
have shown that the magnetic moment is quite small 
(w 0.2/iB/Fe at p = 2.4 GPa). A linear exrtrapolation 
of the moment with pressure would lead to a moment of 
« 0.35/iB a.t p — 4.4 GPa. Therefore, we analyzed the 
neutron data using a theoretical model considering the 
two proposed magnetic structures Ki and Kjj. For both 
structures the magnetic peaks are hidden in the back- 
ground. The simulated diffraction patterns for the stru- 
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FIG. 8: (color online) (a) Difference of the neutron diffraction 
spectra of FeSei-^ taken at T = 5K and 150 K. Only noise 
is seen, except at the position of the nuclear peaks due to 
different Debye- Waller factors. Simulations of the magnetic 
structures (b) my{Ki)-type and (c) my{Kji)-type to the mea- 
sured neutron diffraction spectra at T = 5 K for a moment 
of 0.5 per iron atom. The possible magnetic diffraction 
peaks are hidden in the background signal for all magnetic 
structures proposed. The peak in (b) and (c) indicated with 
* is a temperature independent feature of the pressure cell. 



tures with the magnetic vector Kj and Ku are shown 
in Fig. [Sjj and c. The largest possible magnetic moment, 
that is not seen due to the high background of the pres- 
sure cell is estimated to ~ 0.5 — 0.7 /xb per iron atom (de- 
pendent on the magnetic structure). The simulations of 
the estimated structures are in agreement with the muon 
stopping site calculations that show a very low magnetic 
moment per Fe atom. 



VII. PHASE DIAGRAM 

Figure |9] summarizes the results obtained in this study 
in a phase diagram. At low pressures below p < 0.8 GPa 
the samples are superconducting only and show an in- 
crease of Tc from ~ 8K at ambient pressure to '-^ 13K 
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at ~ 0.8 GPa. At higher pressures static magnetic or- 
der is estabhshed below Tn > Tc that first competes and 
coexists with superconductivity, and at higher pressure 
(p > 1.2 GPa) it only coexists with superconductivity. In 
the intermediate pressure range (0.8 < p < 1.2 PGa) the 
competition is evident from two observations: (i) as a 
function of pressure is suppressed as soon as mangetic 
order appears, leading to the local maximum of Tc at 
p ~ 0.8 GPa. However, the superconducting volume 
fraction remains to be 100%. (ii) the magnetic order, 
that is established above is parially (or even fully)S 
suppressed by the onset of superconductivity. This is 
seen by a decrease of the internal magnetic field i3int(0) 
and a decrease of the magnetic volume fraction when the 
samples enter the superconducting state (see Fig. |4] and 
Ref. 8). For p > 1.2 GPa magnetism is fully established, 
and both Tn and the magnetic moment increase with 
increaing pressure. Interestingly, the onset of magnetic 
order and the simultaneous rapid increase of the Fe mag- 
netic moment coincide with a drastic change of th e Se 
height above the Fe plane that starts also at ~ 1 GPa.'^'^Sl 

The appearance of antiferromagnetic order has also 
been seen by NMR measurements.^ An increase of 1/TTi 
close to Tc is observed at low pressures {p — and 
0.7 GPa) indicating antiferromagnetic modes of spin fluc- 
tuations that are strongly enhanced towards Tc. This 
leads to the conclusion that FeSei_a; is in close proximity 
to a magnetic instability. At higher pressures (at 1.4 GPa 
and 2.2 GPa, i.e. where /iSR observes static magnetic or- 
dering) the 1 /TTi data reveal a broad hump significantly 
above Tc. Furthermore, the integrated intensity of the 
NMR signal begins to decrease at about 34 K at 1.4 GPa 
and at about 50 K at 2.2 GPa, in excellent agreement 
with the /iSR data. The disappearance of the NMR sig- 
nal below a peak of 1 /TTi is a characteristic signal for a 
magnetic phase transition with a (nearly) static magnetic 
hyperfine field with a broad distribution.^ 

Keeping in mind that the superconducting volume 
fraction is ~ 100% for all pressures measured and that the 
magnetic volume fraction reaches ~ 100% at p > 1.2 GPa 
indicates that both ground states coexist in the whole 
sample volume. The data do not show any signature 
for macroscopic phase separation into superconducting 
and magnetic regions larger than a few nanometers, as 
observed e.g. in Bai_2;K^Fe2Asj2ll or LaFeAsOi_a;F2..'2Sl 
No sublattice is present which could order magnetically, 
while the superconducting FeAs layers are not magneti- 
cally ordered, as e.g. observed in Cellll or SmlllLf^iEIl 
These observations point rather to an atomic scale co- 
existence of the order parameters as it is seen e.g. in 
FeTei_a;Sea;'^ or Ba(Fei_j;Co2:)2As2.^^ Furthermore, it 
seems that the two ground states stabilize each other with 
pressure as Tc, Tn, and -Bint(O) are increasing in paral- 
lel with increasing pressure. Comparing FeSci-^: with 
the newly discovered i?Fe2_2;Se2 (245) system in which 
superconductivity and magnetism coexist rises the ques- 
tion, whether magnetic order in FeSei_a; under pressure 
is of similar origin as the one in the 245 system.'^'* In 
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FIG. 9: (color online) Pressure dependence of the super- 
conducting transition temperature Tc, the magnetic ordering 
temperature Tn, and the superconducting and magnetic vol- 
ume fractions of FeSei_a;. The superconducting volume is 
100% for all pressures investigated, determined from ac sus- 
ceptibility and muon spin rotation experiments of FeSei-^;. 
The data obtained in this study are plotted together with the 
data from Refs. [8| and 1111 The Tc and Tn lines are guides 
to the eye and SC, M, and PM denote the superconducting, 
magnetic and nonmagnetic states of the samples, respectively. 



the latter system the superconducting transition temper- 
atures reaches Tc ~ 32 K and superconductivity seems to 
coexists with magnetism occuring at Tn ~ 500 K with a 
rather large magnetic moment of 3/is per Fe atom.l^Sl 

Knowing that FeSei_j: is a two gap superconductoi'^'^ 
a possible scenario of an atomic scale coexistence of 
superconductivity and magnetism has recently been 
proposed by Vorontsov et alW^^ and Gvetkovic and 
Tesanovic.SH They proposed a region in which supercon- 
ductivity and magnetic order can coexist. Here, the mag- 
netic order can be commensurate only in a rather small 
parameter range where the Fermi surface nesting is not 
perfect. The bands are supposed to have an elliptical 
shape, and the chemical potential is supposed to shift. 
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VIII. CONCLUSIONS 

The pressure dependence of the superconducting and 
magnetic properties of FeSei-^; were studied by means of 
ac and dc magnetization, as weU as zero field /iSR tech- 
niques. It is shown that independent on the preparation 
procedure the samples are bulk superconductors up to 
a pressure of p ~ 2.4 GPa. The superconducting transi- 
tion temperature Tc increases with increasing pressure. 
However, the increase is non linear: Tc exhibits a local 
maximum at 0.8 GPa and a local minimum at 1.2 GPa. 
At pressures higher than ~ 0.8 GPa static magnetic or- 
dering occurs below the Neel temperature Tn > T > T^. 
In an intermediate pressure range where Tc is decreasing 
(0.8 < p < 1.2 GPa) the magnetic order is incommen- 
surate and competes with superconductivity.'^ Only at 
p > 1.2 GPa when magnetic order is fully established, the 
magnetic order is commensurate and magnetism occupies 
the full sample volume, coexisting with superconductiv- 
ity on an atomic length scale. Muon stopping site calcula- 
tions reveal only one stopping site of the muons along the 
Se - Se connection and a small pressure dependent mag- 
netic moment with a value of ^ 0.2/iB at p ^ 2.4 GPa is 
found. A recent Mossbauer study reported no magnetic 
order in FeSei _ a: .^^ However, the samples were prepared 
in a slightly different way. Following carefully the prepa- 
ration procedure used in the Mossbauer study and inves- 
tigating these samples by means of /iSR, clear evidence 
of magnetic order in the system is observed,!^ in contrast 
to the Mossbauer results.'^ 

Different magnetic structures based on the muon stop- 
ping site calculations and a symmetry anlaysis are pro- 
posed and tested. The neutron diffraction measurements 
did not reveal any magnetic Bragg reflections because 
the magnetic moment seems to be too small. Thus, only 
speculations about the magnetic structure are possible. 
It is most probably very similar to the magnetic struc- 
ture of the LaFeAsO family of Fe-based superconductors, 
since the FeSei_a; layers resemble the FeAs layers in the 
i?llll system. 

Both superconductivity and magnetism are stabilized 
by pressure. This is evident from the simultaneous in- 
crease of Tc, Tn, and i?int(0) and the related magnetic 
moment /x with increasing pressure. It remains to be 
seen whether this peculiar behavior influences or even 
helps to clarify the pairing mechanism in the Fe-based 
superconductors. 
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A. Magnitude and symmetry of dipole fields from 
Fe subsystems at the muon site 



The symmetry analysis was done assuming that the 
overall distribution of the magnetic fields in the mag- 
netic unit cell has the same symmetry as the magnetic 
order parameter. In order to find the orientation of the 
magnetic field at the muon site, an artifacial magnetic 
moment is ascribed to this site. The corresponding set 
of magnetic degrees of freedom forms the magnetic rep- 
resentation for some positions (Wyckoff positions) . The 
magnetic representation is transferred into an irreducible 
representation after making a standard decomposition. 
After that it is possible to analyze the possible symmetry 
of the mangetic moment (i.e. staggered magnetic fields) 
at the muon site. 

The magnetic order parameters consist of Fourier com- 
ponents of respective magnetic propagation vectors Ki of 
the a sublattice magnetic moments m^\Ki) {a ~ 1, 2): 

k{Ki) =i (m'^\Ki) - mf\Ki)^ ; I = 0,1,11,111. 

(5) 

The nonzero components of respective magnetic mo- 
ments at the muons sites have the form: 

U{Ki) J- (B['\Ki)-Bf\Ki)) ; I = Q, I, II, III. 

(6) 

Here b[°'\Ki) is the i-cartesian component of a mag- 
netic field at the muon site a (a = 1,2) with Ki type 
symmetry. 

In Table [H] the result of the symmetry analysis is pre- 
sented. Here, the enumeration of the irreproducible rep- 
resentations Ti is given in accordance with the Kovalev 
notation.l^ It shows the symmetry of the iron magnetic 
order parameter rjii and U, and the corresponding mag- 
netic fields at the muon sites Mi and Li for the four 
magnetic propagation vectors Ki. Due to the high local 
symmetry of the muon sites some directions of the iron 
magnetic structure cannot create a magnetic field at the 
muon sites. Thus, the observation of /iSR signals (os- 
cillations in the /LtSR time spectra, see Fig. [s]) at high 
pressures in FeSei-a, evidences that the magnetic struc- 
ture has a certain direction and a certain arrangement 
of exchange interactions [i.e. different type of exchange 
order). 
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TABLE II: Symmetry of iron magnetic order parameters (m^ and h that are a linear combination of the sublattice moments, 
see Eq. ([5|), and the corresponding magnetic fields at the muon sites (Mi and Li, see Eq. Q) in FeSei_.a; for the four possible 
propagation vectors of magnetic ordering Ki (l = 0, J, //, ///). The enumeration of the irreproducable representations (IR) 
Ti is given in accordance with the Kovalev notation!^ 
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The analysis of the magnitude and the symmetry of obtained in Eqs. ([7]) and (|8]). There the magnetic fields 

the dipole fields for the possible propagation vectors of are given in mT, and the basis functions (m and I) in the 

magnetic ordering Ki (l — 0, /, //, ///) at the muon units of /^s. 
site of the Fe subsystems in FeSei_a; leads to the results 



For 4 GPa the following results were obtained: 

^B,{Ki)\ /O \ fm^iKi) 

By{Kj) = 354.6 my{Ki) 
^B,{Ki)J \0 354.6 / \m,{Ki) 

^B^{Kii)\ /O \ (m^{Kji) 

ByiKii) = -334.3 my{Ku) 
,B,{Kii) \0 -334.3 / \m,{Kii 




'B,iKo)\ /106.2 \ m^iKo) 

By{Kn) = 111.0 my{Ko 
^B,{Ko)J \ 439.7/ \m,{K„) 

^B.,{Kiii)\ /-217.0 \ (m,{Kni)\ ( 476.1 0\ (UKmf 

By(Kiii) \ ^ \ -222.7 my{Kui) + 476.1 ly{Kjjj) | (7) 

,B,{Kjjj) \ 439.7/ \m,{Kiji) \ 0/ 



For 9 GPa the following results were obtained: 



^B,{Ki) 
ByiKi) 

'B,{Kii 
By{Kn) 
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"yiKlII 
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(8) 
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